Estimation of the large order behavior of the plaquette 
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The universality of vacuum condensate can be exploited to relate the infrared renormalon caused 
large order behaviors of different processes. As an application the normalization constant of the 
large order behavior of the average plaquette is estimated using the Adler function. 



As is well known the perturbative expansion in weak coupling constant in field theory is in general an asymptotic 
expansion, with perturbative coefficients growing factorially at large orders. There are two known sources for this 
behavior. One is the factorial growth of the number of Feynman diagrams at large order, which may be understood 
using the instanton technique [l| . The other is the renormalon in which certain types of Feynman diagrams give rise 
to the large order behavior via their infrared (IR) or ultraviolet (UV) behavior of Feynman integrals (For a review see 
0). These renormalons cause singularities in Borel plane whose properties can be studied by operator insertions for 
UV renormalons and operator product expansion (OPE) for the IR renormalons 0,01 ■ In quantum chromodynamics 
(QCD) the Borel summation of the asymptotic series by IR renormalon is inherently ambiguous, manifested by the 
presence of singularities on the integration contour. This ambiguity in Borel summation is supposed to be cancelled 
by the corresponding ambiguity in the vacuum condensates of the OPE. While this has not been proven there is 
support for it from two-dimensional nonlinear ct- models in solvable large-N limit, where the ambiguity in imaginary 
part of the condensate is correlated with the contour choice of the Borel summation Indeed, the nature of the 

renormalon singularities can be obtained via this cancellation of the ambiguities 0| . The purpose of our paper is to 
use this idea of ambiguity cancellation to relate the large order behaviors of different processes. 

Consider a real quantity G{as) that has an OPE expansion 



where Co denotes perturbative contribution, Oi is the operator for the first power correction, and the suppressed are 
the higher dimensional operators. For simplicity, the dependence on dimensional parameters in the Wilson coefficients 
and condensate are also suppressed. The Borel summation of the perturbative series is ambiguous, which appears as 
contour dependent imaginary term that is to be cancelled by the ambiguity in the condensate {Oi ) . This means that 



where denotes the Borel summed of Co, must be process independent, since the condensate, being a vacuum 
property, should be universal, depending on no particular process. We note that when comparing ^ between two 
quantities the Wilson coefficients are to be computed in the same renormalization scheme, unless the condensate is 
scheme independent. Since the ambiguity is proportional to the normalization constant of large order behavior, this 
implies that the large order behaviors of the quantities that have the OPE ([1]) with common condensate (Oi) are all 
interrelated. To be specific, assume Co has perturbative expansion 



G(a,) = Co(a,)+Cl(a,)(Ol) + ••• 
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This can be expressed in Borel integral as 
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with the Borel transform given by 
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which is expected to have a finite radius of convergence, and /3o is the one loop coefficient of the beta function given 
below in ([5]). The above mentioned cancellation of ambiguities demands the Borel transform have the singularity of 
the form 

G(6) = ^^-^^^(1 + (!?(1-VM), (6) 

where bo is determined by the dimension of the operator and v by the renormalization group equation for the 
condensate and are given as [3, 0] 

n nfSi 71 
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where n is the dimension of Oi and Pi are the coefhcients of the QCD beta function 

2dag{fi) 2 s 
PQCD(as) = fJ- = -Poas ~ Pitts . (8) 

and 71 is the coefficient at O^a^) of the anomalous dimension of Oi. The large order behavior is determined by the 
singularity and is given by 

The singularity causes the Borel integral depend on the choice of the contour, rendering the integral ambiguous. 
Taking the contour along the positive real axis on the upper half plane, the ambiguity, given by the imaginary part 
of the Borel integral 



C^^'i^s)^^ e-'/P°^^G[b)dh, (10) 

Po 



/ 0+is 

where e denotes a positive infinitesimal, is obtained as 

ImC^'^ias) = AAsin(i.^)r(-i.)(6o//3o)'+''e-''«/^«"=a7''(l + 0(a,)) . (11) 

This imaginary part is to be cancelled by that of the condensate, hence 

ImCo^i^(a,) + Ci(a,)Im(Oi) = , (12) 

which means lTaC^^{as) / Ci{as) is process independent. Since the normalization is proportional to the ambiguity 
this allows one to interrelate normalizations among different processes, and also shows that the normalization must 
be proportional to the leading order coefficient of the Wilson coefficient Ci . 

As an application, let us consider the average plaquette and the Adler function. Both have the gluon condensate 

(G^) - -(^^^f^G?.) (13) 

as the leading operator for power correction, hence the large orders of these can be related. The OPE for the average 
plaquette is given by 

P(/3) = (1 - iTrUc) = ^o(a.a) + Z{a,a){G^)a^ + 0{a^) , (14) 

where 

^(/3) = ^(1 + 0(«sd)), (15) 
a is the lattice spacing, and — 3/2tt(3 denotes the bare coupling. The OPE for the Adler function 
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where 

with ~ and a flavor nonsinglet vector (or axial) current, is given by 



n(Q2) = ' / d42.g.9^(o|TJ^(x)J'^(0)|0) , (17) 



D{asm = DoiMQ)) + DiiasiQ))^ + 0(1/Q') (18) 

where 

D,{a,) = ^{l + 0{a,)). (19) 

Since we are interested in QCD with no light quark flavors to compare with the average plaquette of pure Yang- 
Mills theory, we assume that the quarks composing the current are massive so that they do not contribute to IR 
renormalon but still satisfy rn^^^^^ <^ to make the OPE ^TE\\ valid. In this limit, quark bubbles should drop from 
the renormalon diagrams and the only quark lines are those contracting the currents. 

Since n = 4 and 71 = for the gluon condensate (IT^ the renormalon singularity for the plaquette and the Adler 
function can be written, respectively, as 

m « (,_-^°)i+^ (l + g(l-V2)) (20) 



with 



2(3i _ 204 
^"121 



2 - .0. ■ (21) 



Now the ambiguity cancellation between the Borel summed perturbative contribution and the gluon condensate (jl3p , 
along with the renormalization scheme independence of the gluon condensate by the trace anomaly iS] , gives 

Applying the formula (|lip to the Borel integral with the singularities (PU]) we get 



- ^(Q«) 





/3o ois{Q) asu' 



(l + 0(a.)). (23) 



Using the relation between the lattice coupling and the MS coupling at N f = 0, where Nf denotes the number of 
light flavors, 

^ — ^2po HaQ) - AiTh + 0{as) , (24) 



where 



we get 



/3o = ^ , h= 0.46820 (25) 

477 



— t 

A/p = £^ A/-MS ^ 28703 AT^ . (26) 

Note that in obtaining this the higher order corrections in and ([M]) can be safely ignored, since the ratio on the 
left hand side of ([^ is independent of the strong coupling and so they must cancel out. Thus, (pSj) is exact. 
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We now turn to the computation of the normahzation constants. The normahzation constant of a renormalon can be 
computed using the scheme in [HI , which exploits the singularity and analytic property of the Borel transform to 
compute the normalization using the usual perturbative expansion. The result is a convergent series expression of the 
normalization. The speed of convergence of the series depends on the quantity involved as well as the renormalization 
scheme. For instance this yields a rapidly converging series for the static interquark potential or the heav y q uark 
pole mass, rendering the normalization to be evaluated accurately with the first few orders of perturbation [l2l IT3l|. 
Recently, the estimation of the normalization was confirmed by numerical simulation [l^ [Tsj . Numerically, the series 
for the normalization for the plaquette does not converge well at the orders known so far and so it cannot be obtained 
through the scheme. On the other hand, the scheme yields a converging series for the Adler function. 

With the Borel transform (PH)) the normalization Ad is given by 

AAd = R{2) (27) 

where 

R{b)^D{b){l-b/2y+'' (28) 

To express R(2) in a convergent series form the singularity at b — 2 must conformally be mapped so that it becomes 
the nearest singularity to the origin. Since the nearest singularity in 6-plane is the UV renormalon at 6 = — 1 we may 
use a mapping like 

^ = TTb ^''^ 

which maps the singularity at 6 = 2 to one at zq = 2/3, which is the nearest one on z— plane. On z-plane the 
normalization can be written as 

AAd = i?(6(^o)) (30) 

i=0 

where the series is now convergent. The coefficients can be computed from the perturbative series for Dq. 
Dq in MS scheme to five-loop is given as [IB, [13] 

DoiasiQ)) =as + dial + dao^ + ^30^ , (31) 

where as = as{Q)/TT and, at Nj = 0, 

di = 1.98571 , d2 = 18.2427 , = 135.792 . (32) 
The corresponding Borel transform is given as 
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TT 
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(33) 



with which ([30|) gives 



AT^s = - (1 - 0.41393 + 0.08069 + 0.23598) - ^1^2^ (34) 

TT TT 

The series converges well up to four-loop order but jumps at five-loop. This jump is typical of the series for a singular 
function and just may refiect the singular nature of R(b{z)). Note that R(b(z)) is still singular at z = zq, for ^ is a 
fractional number, but being bounded its series is guaranteed to converge at zq; Nevertheless, the convergence can be 
bumpy, unlike for the series of a smooth function. It is also interesting to see the behavior of the normalization at 
differing Nf. For the first few nonzero flavors we have 

r (1 -0.40840 -(-0.01607 -f0.18119)/7r for iV/ = 1 , 

A/]^^ = S (1 - 0.39613 - 0.06313 -1- 0. 11121)/7r for Nf^2, (35) 
[ (1 -0.37421 - 0.16118 -h0.01910)/7r for A^/ = 3 , 

which shows a better convergence at increasing flavor numbers, a behavior that was already observed with Adler 
function of electromagnetic current [loj . 
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FIG. 1: Scale dependence of the normalization constant at jJ. = Q of the Adler function. 

Now, taking the five-loop contribution as the uncertainty in the estimate we conclude 

_^^^a90±a24^ (36) 

TT 

from which we get the normalization for the average plaquette: 

25833 ±6889 

Afp = . (37) 

TT 

Considering the jump at five loop the uncertainty in this estimate can be too optimistic. To avoid such underestimation 
of the uncertainty we may also look at the renormalization scale dependence of the normalization constant. Being 
proportional to the gluon condensate the ratio 

is scale independent, which means that the normalization 7V^^ (/x) of the series in powers of as (^) of the Adler function 

scales as 

A/r(A.)=A/-r(Q)(/./g)^ (39) 

Thus the test of scale independence of 

N-^{Q)^M^{ii){Q/f,r (40) 

can give a hint of the reliability of the estimate (1551) . In Fig. [T] we see that Q is close to the scale of minimal 
dependence, and considering the variation of the normalization about /i = Q it appears the error estimate in p6|) is 
not unreasonable. 

At this point it may be appropriate to isolate the universal portion of the normalization constants that is process- 
independent, by writing the normalization as 

Af^ = coA/gF (41) 

where cq denotes the leading order coefficient of the process-depedent Wilson coefficient for the operator . From 
the estimate of the normalization for the Adler function and (IT9l) we obtain 

/■MS 



1.35 ±0.36 



^MS ^ — . (42) 

The process-dependence of the normalization comes via Co, a short-distance quantity, and Af^^ is the process- 
independent part of the normalization, which may be regarded as the long-distance contribution and an intrinsic 
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FIG. 2: Renormalon (solid) vs power law (dashed) behavior. 



property of the renormalon like the strength v or the position of the renormalon singularity. That the process- 
dependence comes only via a short-distance quantity should not be surprising, considering that in IR renormalon 
diagrams the large order behavior arises from bubble chains of arbitrarily long, far-infrared region; hence all process- 
dependence should be a short-distance effect. 

Now note that the normalization constant p7p is for the expansion in ctsu- Fo^' the usual power expansion in 



(43) 



the large order behavior is then given by 



87r7Vpr(i -I- zy) / 33 



(1 + 0(1A)) 



(44) 



The plaquette coefficients were computed in numerical stochastic perturbation theory up to 20- loop orders |18l - l20| . 
At these orders the coefficients grow much faster than a renormalon behavior would suggest, and rather follow a power 
law. The plot (Figl^]) of the renormalon behavior (pi)) and power law [2lj shows they meet at order i ~ 42. This 
may suggest the renormalon behavior would set in at orders around i ~ 40. Recently, the renormalon behavior in 
heavy-quark pole mass was confirmed in numerical simulation of the coefficients to order o?^ [13, [iBl • Our estimate 
of the large order behavior of the plaquette suggests a numerical evidence of renormalon in plaquette would require 
much higher order computations. 
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